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Abstract 

The so-called min-sum algorithm has been applied for decoding lattices con- 
structed by Construction D . We generalize this iterative decoding algorithm 
to decode lattices constructed by Construction D. An upper bound on the 
decoding complexity per iteration, in terms of coding gain, label group sizes 
of the lattice and other factors is derived. We show that iterative decoding 
of LDGM lattices has a reasonably low complexity such that lattices with 
dimensions of a few thousands can be easily decoded. 
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1. Introduction 

Both of the integer programming method and the trellis approach, as 
two main methods for lattice decoding, are impractical in higher dimensions 
[4, 5]. The min-sum algorithm, as an iterative decoding approach, can be 
used in decoding high dimensional lattices. Tanner generalized Gallager's 
low-density parity-check (LDPC) codes to other class of codes defined by 
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general bipartite graphs, called Tanner graphs [7]. For any linear binary 
block code, this construction is based on a parity check matrix of the code. 
Tanner graph construction is used to find a graphical representation of the 
lattice in terms of its linear constraints. The decoding complexity of the 
generalized min-sum algorithm depends on the Tanner graph structure and 
the label code of the lattice. Sadeghi et al. introduced a generalization of 
min-sum decoding algorithm for lattices constructed by Construction D [6]. 
In this work, we will propose another generalization of min-sum algorithm to 
decode lattices constructed by Construction D. Therefore properly selected 
lattices, such as those based on low-density generator matrix (LDGM) codes, 
can be decoded efficiently. The paper begins in the next section with a brief 
discussion about lattice. Section three introduces the generalized version of 
min-sum algorithm to decode lattices constructed by Construction D. The 
decoding complexity and it's bounds for the new generalization of min-sum 
algorithm are discussed in the forth section. The final section is dedicated 
to the paper's conclusions. 

2. Preliminaries 

Let M m be the m-dimensional real vector space with the standard product 
(., .) and Euclidean norm || x ||= (x, x} 1 / 2 . A lattice A is a discrete addi- 
tive subgroup of K m . An n-dimensional lattice is generated by the integer 
combinations of a set of n linearly independent vectors [3]. Any subgroup 
of a lattice A is called sublattice of A and a lattice is called orthogonal if 
it has a basis with mutually orthogonal vectors. The set A* of all vectors 
in the real span of A (span(A)), whose the standard inner product with all 
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elements of A has an integer value, is an n-dimensional lattice called the dual 
of A. Let us assume that an n-dimensional lattice A has an n- dimensional 
orthogonal sublattice A'. If A' has a set of basis vectors along the orthog- 
onal subspaces S = {W / i }™ =1 , the projection onto the vector space Wi de- 
fined as Py/ l and the cross section A Wi defined as A w . — A D W». The label 
group of the lattice is defined as Gi = Pw i /Aw i , which is used to label the 
cosets of A' in A. For any lattice-word x, the label sequence is defined as 
g{x) = (pi(x), . . . ,# n (x)), where &(x) = P Wi + A Wt . The set of all possi- 
ble label sequences, L = g(A) = {g(x) : x £ A}, is the label code of A. 
This set is an Abelian group block code over the lattice-alphabet sequence 
space, G = Gi x ... x G n . Let \Gi\ = g,i and v$ be the generator vector of 
Ay/., i.e., Ayy. = Zv.j. Each element of Gi can be rewritten in the form of 
+ j det(P Wi )Vi/\vi\(j = 1, . . . ,gi - 1). Then the map 

A Wi +jdet(P Wi )p- — > j (1) 

I i I 

is an isomorphism between Gi and Z 5i [6], thus every element of the label 
group Gi can be written as (Z + <2j)vj, where a, = j det(Pwi) / 'det(Awi) ■ 

There exist another efficient method for lattice representation introduced 
by Tanner [4]. Lattices constructed by Construction D have a square gen- 
erator matrix. If B is a generator matrix for A, then B* = (B _1 ) ir is a 
generator matrix for A* (parity-check matrix of A) [2]. This can be applied 
to construct the Tanner graph for the lattice [3]. To construct the Tanner 
graph for a lattice, The Tanner graph construction of linear codes is applied 
to the parity check matrix of the lattice. If Sj denotes the ith column and chj 
denotes the jth row of B* respectively, then the Tanner graph of the lattice 
has the edge (sj, chj) if and only if the lattice- word component (symbol node) 
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Si is contained in (or checked by) the parity-check sum (check node) chj. As 
a result x G Z n belongs to A if and only if B*x T G Z. 



Example 2.1 Consider the following B and B* as the generator matrix and 
the parity check matrix of a 7-dimensional lattice constructed by Construc- 
tion D. 
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The corresponding Tanner graph for this example is illustrated in Fig. 1. 
The black circles denote lattice-word components and the white rectangles 
represent parity-check sums. 



3. Generalized min-sum algorithm 

In this section, we explain how generalized min-sum algorithm for lat- 
tices constructed by Construction Dl can be changed for decoding lattices 
constructed by Construction D. 

3.1. Lattice decoding using min-sum algorithm 

Given a vector y G M n , the lattice decoding problem is to find a lattice 
vector x, such that || y — x || is minimized. Let y = Xir=i^ v *' wnere 




Figure 1: The corresponding Tanner graph for the 7-dimensional lattice discussed in ex- 
ample 2.1. The black circles and the white rectangles represent symbol nodes and check 
nodes respectively. Si denotes the ith column and chj denotes the jth row of parity-check 
matrix B*. 

Vi = (y> v «) / ( v i? v «)- The goal of the decoding algorithm is to search for the 
closest lattice point to the received-word. By definition: 



.r 



(0 



a, + \yi -a,j\, i = 1, . . . , n and j = 0, . . . , \Gi\ - 1 (2) 

where \u\ is the closest integer to u, aj = j det(P Wt ) / det(A w .) and 
denotes the closest point in jth coset of Gi at the ith coordinate of the 
lattice-alphabet sequence space. Then XjVi is the closest vector of Gi to 
j/iVj. The set xW := : j = 0, ...,gi — 1} has gi candidates of the 

ith coordinate of the lattice-alphabet sequence space for every component of 
the received-word. Define the weight as the squared distance between the 
elements of (Z + %)vj and yiVi. 



I (*) - \2 II ||2 

\A - Vi) II Vj || 



(3) 



Considering the alphabet sequence space Z 9l x . . . x Z 9n , we rename j to 
Q where c« G Z ffi . The weight of any valid codeword c = (ci, 



,c n ) Gi 



correspond to lattice-word x = (xci, . . . , Xcn) G A, is defined as 

n 

^y( c ) = ^w w (ci). (4) 

i=i 

Now the problem is to find the minimum weight, min (w y (c)). 

3.2. Min- Sum algorithm for Construction D lattices 

This algorithm includes initialization, hard decision, symbol node opera- 
tion and check node operation. 

Lemma 3.1: For any received vector y = . . . , y n ) G M. n we have 



= det ( A 



(t 



Vi 3 



y% 3 



.det(k Wi ) gj det(A Wi ) g 
where % — 1, . . . , n and j — 0, . . . , gi — 1. 



(5) 



Proof. Proof is given in [6]. □ 

1) Initialization: Let y = Ym=\ = Y^i=i Vi e i e ^ e the received 
vector. An initial weight is assigned to each node: 



i = (u y M,---,Uyi(9i ~ 1)) (6) 



2) Iteration: In the iteration step all weights are alternatively updated 
to find the lattice-word. 

2.1) Symbol node operation: The intermediate weight, which de- 
notes the symbol-to-check outgoing weight, computed as follows: 

u yi ,ch(k) := u y .(k) + ^2 u ch >, yi (k), < k < gi - 1 (7) 

ch'eQ 
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2.2) Check node operation: The intermediate weight, which de- 
notes the check-to-symbol outgoing weight, computed as follows: 

Uch, Vl ( k ) '■= min 5^ u y c h (k'), Q<k'<g v -1 (8) 

_ (i) y^ech 

where Q denotes the set of check equations and Q c h denotes the 
set of all valid configurations that satisfy the chth check equation. 



3) Termination: For every symbol node, all incoming messages add to 
its initial weight to obtain final weight as follows: 

Fu y .(k) := u yi (k) + ^ u ch i >yi {k). (9) 

ch'eQ 
ViSch' 



The goal of the min-sum algorithm is to find a vector x = (x Cl , • • • , x c „ j 



which x^' G {xq j , . . . , a^-i} and the index k, obtained as follows: 



k = arg(unn Fu Vi {k)). (10) 

0<k<gi — i 

The iteration will stop when the selected vector x is a lattice- word, i.e., 
x satisfies all parity check equations or reaches the maximum number of 
iteration. 

4. Decoding complexity 

In each Termination step for every lattice- word component (symbol node) 
with label group size gi, the number of comparisons is (ft — 1. Therefore the 
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total number of comparisons is: 

n 
i=l 

In each Iteration step, the number of operations for each symbol node and 
check node computed separately: 

a) For each symbol node and for each edge, the number of summation is 
gi(d y . — 1), where d Vi denotes the number of edges for each symbol node. 
Thus the total number of summation would be: 

n 

^Mfctefc-l). (12) 

1=1 

b) For each check node, at most, x ... x g id ^ comparisons should be 
made. For each outgoing message, u c h uyk (j), and for each j, the number of 
summation is d c } H — 1, where d c ^ denotes the number of edges of each check 
node. Since there are d c h t edges, the number of summation will not exceed 
dchXdch, ~ l)(<?ii x • • • x 9i d h ) summations for each check node. Then the 
total number of operations in all check nodes in each iteration is at most: 

n 

^2d chi (d cht - l)(g h x ... x g idch ). (13) 
i=i 

Eqs. (11), (12) and (13) show the dependency of the decoding complexity on 
the size of label groups. In each iteration the total number of operations is: 

n 

[J2 9i + 9id yi (dy t -i) + d chi {d chl - l)(g h x . . . x g idch )) - n. (14) 

i=l 

The next Corollary follows from counting the number of operations and in- 

1 /2 

equality g^ > (7(A)7(A*)) , where 7(A), denotes the coding gain of the 
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lattice [1]. 



Corollary J^.l: (Bounds on decoding complexity) Let A* be the dual of A 

1 /2 

and 7 = ( 7 (A) 7 (A*)) 1 . Also assume that A has a Tanner graph with n 
symbol nodes and n check nodes for which g,i < g, d y < d y . < d^ ax and 
d c h < d c hi < d™£ x {i = 1, . . . ,n). The upper bound of decoding complexity 
per iteration is: 

n(gd% a *{d% a * - 1) + g d ™ x d2 ax (d2 ax - 1) + g - l) (15) 

and the lower bound per iteration is: 

n(jd y (d y - 1) +Y ch d ch {d ch - 1) +7- l). (16) 

The proof is a direct consequence of Eq. (14) and the fact that gi > 7 which 
has been shown in [1]. 

This corollary shows that the decoding complexity per iteration, grows lin- 
early with the lattice dimension, n, but has power law dependence on the 
check nodes degree. 

5. Conclusion 

In this work the min-sum algorithm is generalized to decode lattices con- 
structed by Construction D. It is shown that the upper and lower bounds 
of decoding complexity depends on lattice parameters like label group sizes, 
coding gain, check nodes and symbol nodes degree of Tanner graph. It is also 
shown that the decoding complexity grows linearly with the lattice dimen- 
sion, n, but has the power law dependence on the check nodes degree. The 
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analysis of decoding complexity confirms the usefulness of LDGM codes for 
the lattice construction. It is worth mentioning that the presented decoding 
algorithm can be used to decode other constructions of lattices from linear 
codes. 
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